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We review standard subtraction, as a method to compute cross sections at NNLO accuracy. 
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1. Introduction 

High-energy physics will enter a new era of dis- 
covery with the start of LHC operations in 2007. 
The LHC is a proton-proton collider that will 
function at the highest energy ever attained in 
the laboratory, and will probe a new realm of 
high-energy physics. The use of a high-energy 
hadron collider as a research tool makes substan- 
tial demands upon the theoretical understanding 
and predictive power of QCD, the theory of the 
strong interactions within the Standard Model. 

At high any production rate can be ex- 
pressed as a series expansion in as- Because QCD 
is asymptotically free, the simplest approxima- 
tion is to evaluate any series expansion to lead- 
ing order in as- However, for most processes a 
leading-order evaluation yields unreliable predic- 
tions. The next simplest approximation is the 
inclusion of radiative corrections at the NLO ac- 
curacy, which usually warrants a satisfying as- 
sessment of the production rates. In the previ- 
ous decade, a lot of effort was devoted to de- 
vise process-independent methods and compute 
rates to NLO accuracy. In some cases, though, 
the NLO corrections may not be accurate enough. 
Specimen cases are: the extraction of as from fit- 
ting the predictions to the data, where in order to 
avoid that the main source of uncertainty be due 
to the NLO evaluation of some production rates, 
like the event shapes of jet production in e^e~ 
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collisions, only observables evaluated to next-to- 
next-to-leading order (NNLO) accuracy are con- 
sidered P ; open 6-quark production at the Teva- 
tron, where the NLO uncertainty bands are too 
large to test the theory |21 vs- the data Higgs 
production from gluon fusion in hadron collisions, 
where it is known that the NLO corrections are 
large gEj, while the NNLO corrections P7I5] . 
which have been evaluated in the large-mt limit, 
display a modest increase, of the order of less 
than 20%, with respect to the NLO evaluation; 
Drell-Yan productions of W and Z vector bosons 
at LHC, which can be used as "standard can- 
dles" to measure the parton luminosity at the 
LHC 

In the last few years the NNLO corrections 
have been computed to the total cross sec- 
tion |6I13| and to the rapidity distribution |14l 
[T5] of Drell-Yan production, to the total cross 
section for the production of a scalar |6l7l8j 
and a pseudoscalar |16I17| Higgs from gluon fu- 
sion as well as to a fully differential cross sec- 
tion |18I19| . and to jet production in e+e~ col- 
lisions [20121122] . The methods which have been 
used are disparate: analytic integration, which 
is the first method to have been used jTS], can- 
cels the divergences analytically, and is flexible 
enough to include a limited class of acceptance 
cuts by modelling cuts as propagators |7I14I15[ 
ITT) : sector decomposition |18I20I23I24I25I26I27| . 
which is flexible enough to include any accep- 
tance cuts and for which the cancellation of the 
divergences is performed numerically; subtrac- 

tion li>:jl:jHli>kl:ji)l:^0l;^ll:^:jl;^:^l:^4l . for which the 
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cancellation of the divergences is organised in a 
process-independent way by exploiting the uni- 
versal structure of the infrared divergences of a 
gauge theory, in particular the universal struc- 
ture of the three-parton tree-level splitting func- 
tions |35l36l37l38l39l40j and the two-parton one- 
loop splitting functions j41l42l43l44l45j . 

The standard approach of subtraction to 
NNLO relies on defining approximate cross sec- 
tions which match the singular behaviour of the 
QCD cross sections in all the relevant unresolved 
limits. For processes without coloured partons 
in the initial state, we constructed in '34' sub- 
traction terms which regularise the kinematical 
singularities of the squared matrix element in all 
singly- and doubly-unresolved parts of the phase 
space. Thus, the regularised squared matrix ele- 
ment is integrable over all the phase space regions 
where at most two partons become unresolved. 

2. Subtraction scheme at NNLO 

The NNLO correction to any m-jet cross sec- 
tion reads 



,NNLO 



rn+2 

J^VV T 



Jm+1 

(1) 



The three integrals in Eq. are separately diver- 
gent, but their sum is finite for infrared-safe ob- 
servables. As explained in [331, rewrite Eq. ^ 
as 



,NNLO 



.^NNLO 
dCTm+2 

m+2 



/ 1 NNLO 



(2) 



where each integral in Eq. ^ is finite by construc- 
tion. Here we will focus on the subtractions that 
regularise doubly real emission, so we recall only 
that 



1 NNLO 

dcr,„+2 



+2 ^"1+2 



dR.R..AQ T 
^m+2 



(3) 



where dcr^'^'^^ and dcr^'^l'^^ regularise the 
doubly- and singly-unresolved limits of da^^2 re- 
spectively and Acr^2^^ accounts for their over- 
lap. 

3. Subtractions for doubly-real emission 

The cross section Aa^^2^ the integral of the 
tree-level squared matrix element for (m -1- 2)- 
parton production over the (m -I- 2)-parton phase 
space 

da^^,=d0("+^)|A^|°V2p- (4) 
The counterterms may be written symbolically as 

I RR.A, 

d^™+2 



[d|,(i)]A|A^i°V2^ (6) 



and 



m+2 1 • (''') 



In this contribution we define explicitly the 
-^i|-^m+2p term, the terms A2\M^°n+2? and 
^i2|A^m+2p will be presented elsewhere. 

The singly-singular subtraction ^i|A^|^^2p is 



MMt\2? =E (^^- - ^^r^r) +T.^r. (8) 



The singly-collinear term in Eq. ((SJ reads 



(9) 



where Pf'^f = ^ff {^i,r, ^r,i, k±^ir', e) is the 
Altarelli-Parisi splitting function. The momen- 
tum fractions Zi^r and Zr^i are defined as 

Zi.r — — and Zr^i — — , (10) 

■S(ir)Q ' ■S(ir)Q 

and the transverse momentum k±^ir is given by 

^l,ir = Oi,rPr " Qr,iPi + {Zr,i - h,r)pt , (U) 

with 



3(0) 



6^0 



y.irS{ir)Q 



(12) 
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We used the abbreviations s^q = 2pi ■ Q, s^q = 
2pr ■ Q and S(ir)Q = s^q + s^q above. 

The m + 1 momenta entering the matrix ele- 
ments on the right hand side of Eq. (O are 



I - air 

In Eq. H13() n ^ i,r and 



1 - a,> 



2s 



(13) 



(14) 



with Q^^ the total four-momentum of the incom- 
ing electron and positron and s — Q^. 
The singly-soft term is 

if r is a gluon, and iS^ if r is a quark or anti- 
quark. The m + 1 momenta entering the matrix 
element on the right hand side of Eq. p5|l are 

P>;, = A^AQAQ-Pr)/>^rMJXr), n^r, (16) 



where = ^ 1 — s^q / s and 

2{K + KY{K + K)^ iKf'k^, 



{K + Kf 



(17) 



The matrix [K, K] generates a (proper) 
Lorentz transformation, provided — ^ 0. 
In Eq. (|15|l Sik{r) denotes the eikonal factor 



Sik{r) 



2s 



ik 



(18) 



and the sum in Eq. H15|l runs over the external 
partons of the (m -|- l)-parton matrix element. 
The soft-collinear subtraction is defined by 

Sir ^r,i 

if r is a gluon, and CirSr = if r is a quark or 
antiquark. The momentum fractions are given 
by Eq. (|10|l . The correct variables in the matrix 
element in the soft-collinear limit are those that 



appear in the soft limit . Thus the m + 1 mo- 
menta entering the matrix elements on the right 
hand side are given by Eq. pt)|l . 

The momentum mappings introduced in 
Eqs. H13|l and (|16|l in addition to conserving total 
four-momentum, both lead to exact phase-space 
factorisation in the form 



,(m+2) _ J^(m+1) rj„(l) 



(20) 



where the m+1 momenta in d0'^'"+^^ are precisely 
those defined in Eq. ifT^ or Eq. |[TB|) . 
The explicit expressions for [dp*^^^] read 



[dp 



(1)1 = 



J 



d'^p, 



and the Jacobian factors are 

(1 - a,rr^^-''>-' 0(1 - <^^r) S^^q 



J 



(s ~ -f — Sro)'^ + 4s ~ (s 

V irQ ' V / rtr ^ 



(21) 



:(22) 



(23) 



for the coUinear and soft phase-space factorisa- 
tions of Eqs. H13|l and respectively. 

In Eq. air is understood to be expressed 
in terms of the variable pir , 



a,: 



i^rir + ^irO ^tqY + "^^riri^ ^iro) 



2{s - s^^ (24) 



The analytical integration of the counterterms 
over the factorised one-parton phase-space [dp^^^J 
is now possible. Starting with the coUinear sub- 
traction, notice that k±,ir as defined by Eq. pi|l 
is orthogonal to pir, therefore, the spin correla- 
tions generally present in Eq. @ vanish after az- 
iniuthal integration Thus when evaluating 

the integral of the subtraction term Cir over the 
factorised phase space [dp^^^], we may replace the 
Altarelli-Parisi splitting functions pj"j. by their 

azimuthally averaged counterparts Pf'\ , so 



;dp(i)] C^r = a^C„{y-ro'^rn, e) T}r \M 



(0) ,2 
m+1 1 I 



(25) 
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where y~g = 2p„ ■ Q/Q'^ 
and 



(26) 



2e 



Sir T 



(27) 



The function Gir{y^^Q]m, e) depends on the 
momentum of the parent parton and the flavours 
of the daughter partons. ExpHcitly 



'2e 



+ (l-e)lW(x;e) 



C,^{x-m,e)^^x-^'[l[^J{x-e) 



(28) 



1-e 



(lW(x;£)-e(x;£))], (29) 



and 



Ggg {x- m, e) = 2x-^' [ 2 (l^^^ (x; e) - l^^^ (x; e)) 
+ lW(x;£)-l(^'(a;;e) 



(30) 



The analytical formulae for the Jm ■* (a;; e) func- 
tions involve Beta functions, the ^Fi hypergeo- 
metric function as well as the first Appell func- 
tion Fi and a certain generalisation of the last, 
i^P', see Tabled The expansion of Im\x;e) in 
powers of e is performed using the techniques of 
I47I48I to obtain 



(32) 



- 1 


3 

+ Ye- 


2 

— ln(x) 

e 


Cf 


r 2 


-0(£°) 


- 2 


11 

3e 


4 

— ln(a;) 

£ 



The finite parts, not shown here, depend on m 
and can be easily found for any given m using 
the program of Ref . @H| ■ 

Integrating the soft subtraction, the color cor- 
relations of Eq. (|15|) are still present in the inte- 



grated expression 



x|-^!n'M;(..)r (34) 

where 

-87ra,/i''y"[V]^5,fc(r). (35) 

For S,fc {y~f^ , yjg , yj^^ ;m,e) we find 

m(l -£)(!- 26) 

s.fc(ys, yiQ,ykQ;m, s) = 

X B(l - e, 1 - e)B{l - 2e, m(l - e)) 
XK2Fi(l,l,l-e, 1-k), (36) 

where n = yikKyiQVkQ) ^^^^ ^^'^ hypergeometric 
function can be expanded using 

K 2i^i(l, 1, 1 - e, 1 - k) = 



(37) 



The integral of the soft-coUinear counterterm, 
Eq. (Unjl, reads 



with 



1 2£- 



(39) 



J Sir ^i.r 

Performing the integrations, we obtain 



xS(l-e, l-e)B(l-2e,m(l-e)). (40) 
In order to show that the integrated subtrac- 
tion term Ai|A^^^|_2p ^las the correct pole 
structure j46| . we exploit colour conservation 
in Eqs. ijSnil and Indeed, using = 

— X^fc^^i T.iTk, we can combine the soft and soft- 
coUinear subtractions in Eqs. H34|l and (|38|l . 



(41) 



i^r 



{i,k)\ 



r i,k 
i^k 
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Table 1 

Definitions of tlie Im\x; e) functions. {a)m = r(TO + a)/T(a) denotes the Poclihammer symbol. 



/i"'^(x;£) 



B(-2e, l-e)B(-e,l-£) 

X (-2£, -e, 2m(l - e) - 1 - 2£, -2(m - 1)(1 -£);!- 3£, 1 - 2£; 1 - a;, 1) 
+B{1 - 2£, 1 - £)B(1 -£,!-£) 

(1 - 2£, 1 £, 2to(1 - £) - 1 - 2£, -2(m - 1)(1 - £); 2 - 3£, 2 - 2£; 1 - x, 1) 



B{1 - £, 1 - £)B(-£, 2m(l - £))2Fi(2m(l - £) - 1 - 2£, -£; 2to(1 - e) - e;l - x) 



lln\x-e) 



^ ( . j B{k - i + 1 - £, 1 - e)B{i - £, 2m(l ~e) + k-i) 

Fi{i - £, 2m(l - £) - 1 - 2£, k; 2m(l - £) + A; - £; 1 - -1) , fc = 1, 2 



OO CX2 



mi —0 ms— 



mil 7712! 



Expanding (S^/e + Ci^Sr) in powers of e, wc find 

S.fc + C„S, = iln^^ + 0(£°). (42) 
^ yiqVkQ 

We can combine these contributions with the 
collinear functions in Eqs. 131|I - H33|) after using 
colour conservation in Eq. (|25() and find that the 
Iny-.g/e terms in Eqs. and Eq. ex- 

actly cancel and the known structure of the one- 
loop squared matrix elements is reproduced. 

4. Conclusion 

We outlined a subtraction scheme for com- 
puting cross sections at NNLO accuracy using 
the known singly- and doubly-singular limits of 
squared matrix elements. A way of disentan- 
gling these overlapping limits was presented in 
[34^ . Here we discussed how to define the singly- 
singular subtraction terms. We presented exact 
phase space factorisations that allow us the inte- 
gration of the singular factors and the results of 
these integrations. We have also coded Eq. Q for 
the case when d(T^^2 is the fully differential cross 
section for the process e'^e^ qqggg {m = 3) 
and J5 defines the C-parameter. We found that 
the integral of da^j^^ is finite and integrable 
in four dimensions using standard Monte Carlo 
methods. 
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